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Abstract 

Using a merged beams apparatus we have measured the associative detachment (AD) reaction of 
H + H — > H2 + e~ for relative collision energies up to E r < 4.83 eV. These data extend above the 
1 eV limit of our earlier results. We have also updated our previous theoretical work to account 
for AD via the repulsive 2 £g potential energy surface and for the effects at E r > 0.76 eV on 
the experimental results due to the formation of H2 resonances lying above the H + H separated 
atoms limit. Merging both experimental data sets, our results are in good agreement with our 
new theoretical calculations and confirm the prediction that this reaction essentially turns off for 
E T > 2 eV. Similar behavior has been predicted for the formation of protonium from collisions of 
antiprotons and hydrogen atoms. 
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I. INTRODUCTION 



One of the simplest molecular formation reactions is associative detachment (AD) via 



H- + H ->• H 2 + e 



(1) 



This reaction is of interest for fundamental atomic and molecular physics and also because 
it plays an important role in protogalactic and first star formation in the early universe 
jj-yj . Two groups have recently reported measurements of this reaction. Martinez et al. 4J 
measured the thermal rate coefficient at 300 K using a flowing afterglow technique. Our 
group has measured this reaction over a collision energy range from 4 meV to 1 eV using a 
merged-beams method {jUsl,!^]. Our results lie 2.2±0.9 times above those of 4|. The quoted 
uncertainty represents the quadrature sum of the estimated total experimental la confidence 
level for each measurement; and we have also taken into account minor corrections to our 
earlier data which are described below. 

In jo] we hypothesized that this discrepancy is due to an error in the measured rate 
coefficient of 7] for 

H + CT -> HCl + e~ (2) 

which {4] used to determine their neutral H number density and thereby normalize their 
results. Our apparatus is not configured to study reaction (j2j) and test this hypothesis, but 
we have been able to extend our measurements of reaction (Q]) to higher energies and thereby 
provide additional benchmarks for theory. We have also investigated and ruled out several 
possible sources of systematic errors in our previous experimental results. Additionally we 
have more carefully considered the pressure dependence of our detection method. Lastly, we 
have updated our previous theoretical results of {3, 8] to account for AD via the repulsive 
2 E+ state and for the effects at E T > 0.76 eV on the experimental results due to the 
formation of long-lived H2 resonances lying above the H + H separated atoms limit. 

The rest of the paper is organized as follows: Section [Til describes the experimental method 
and the various modifications performed for this work. Section [LTTI discusses the experimental 
uncertainties. Our new theoretical calculations are briefly described in Sec. [TV] In Sec. [V] 
we present our results and compare them to theory. A discussion of our results is given in 
Sec. I VI I and a short summary in Sec. IVHI 
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II. EXPERIMENT 



Here we briefly describe the experiment and the changes relevant to our new results. 
Further details about the apparatus and experimental method can be found in [j], [5], [fl . 

A. Method 

We begin by extracting H~ from a duoplasmatron source and forming a beam with an 
energy of E H - = —e(U s + Uf/2). Here e is the unit charge, U s ~ —10 kV is the nominal 
source voltage, and Uf is a small correction voltage defined below. Using standard ion 
optical elements, we shape, steer, and direct the beam into a photo detachment chamber 
which houses a floating cell biased to a potential Uf. The anion energy inside the floating 
cell is .Eh- — — e(U B + Uf/2) + eUf. Near the center of the floating cell, we cross the anions 
with an infrared laser and convert a portion of the H~ beam into a beam of ground state 
H atoms of energy En = —e(U s — Uf/2). The resulting merged beams exit the floating 
cell, whereupon the H~ beam returns to its initial energy while the H beam energy remains 
unchanged. The beam-beam interaction energy is controlled by varying Uf. 

Shortly after leaving the photodetachment chamber, the two beams enter an interaction 
region of length L. Two beam profile monitors (BPMs) are used to determine the beam- 
beam overlap (£l(z)) within the interaction region, where the z axis is defined by the bulk 
velocity vectors of the co-propagating beams. We also use the BPMs to verify the alignment 
of the beam axes. The relative energy E T between the beams depends, in part, on this 
alignment and is given by 9| as 



£ r = /W — + — W^^ostf • (3) 
\m H - m H V m H -m R J 

Here /x = nin-mn/ (inn- + mn) is the reduced mass of the colliding system; mjj- and me are 
the masses of the H~ and H, respectively; and 9 is the angle of intersection. E T is controlled 
by varying Uf. This merged beams approach allows us to reach collision energies on the 
order of a few meV, limited only by the alignment of the beams, the spread in collision 
angles between the two beams, and the energy spread of each beam. We used geometrical 
simulations 0, 0, 0] to determine the average collision energy (E T ) versus Uf, taking into 
account the spreads in beam energies and angles. 
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Both beams are chopped out of phase in order to extract the signal H 2 generated in the 
interaction region from various backgrounds. Any H2 formed in the interaction region has 
an energy of E^ 2 = E u - + Eu = —2eU s = 20 keV, neglecting the < 3.7 eV kinetic energy 
of the detached electron. At the end of this region, an electrostatic quadrupole deflector 
is used to direct the H~ into a Faraday cup where the current /h- is read and recorded. 
The parent H and daughter H2 beams continue on into a gas cell kept at a helium pressure 
of 2 x 10~ 4 Torr for most measurements. Inside the cell a fraction of the H2 is ionized by 
the stripping collisions forming 20 keV H^. Additionally, stripping of the H beam and 
dissociative ionization of the H 2 can produce ~ 10 keV H + . 

After the gas cell, the neutrals and resulting ions enter the analyzer region involving 



two double-focusing, electrostatic cylindrical deflectors in series [10j and a channel electron 
multiplier (CEM). A hole in the outer plate of the first or lower cylindrical deflector (LCD) 
allows neutrals to pass through and travel into a neutral detector. The neutral particle 
current In, as measured in amperes, is monitored by measuring the secondary negative 
particle emission from the target inside the neutral detector. The voltages on the LCD and 
upper cylindrical deflector (UCD) are selected to transmit the 20 keV signal ions into 
the CEM while rejecting any of the 10 keV H + formed in the gas cell. 

We study reaction ([1]) from the number of ions detected in the CEM. Experimentally, 
we measure the cross section a ad times the relative velocity v r between the H~ and H beams 



convolved with the velocity spread of the experiment. This gives the rate coefficient 



6] 



1 S e 2 vn-vn 

0- st iVHe T^TgT] I U -I K L{il(z)} 

The left hand side average is over the experimental energy spread. On the right side, a st is the 
stripping cross section for H2 on He forming ; Nn e is the gas cell helium column density; 
S is the background subtracted, pressure corrected signal; T a is the transmittance of the 
combined LCD-UCD analyzer; T g is the transmittance of the grid in front of the CEM; rj is 
the CEM efficiency; and v h- and t>H are the velocities of the H~ and H beams, respectively. 

B. Modifications 

The present work uses a current meter with a fast response time which enables us to 
directly measure the H~ current at each phase in the chopping pattern, which is on the 
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millisecond scale, and monitor it throughout each data run. Thus we are able to measure 
the anion current when the laser is on, I^L , over the course of a data run. This is used for 
Jjj- in Eq. (111) . We were also able to monitor the anion current with the laser off, I^L, and 
determine the attenuation factor 

ron 

/ = 1-^. ( 5 ) 

n 

which is needed to extract the background corrected S [6\. This situation is to be contrasted 
with our previous results where, due to equipment limitations, the H current was 

averaged over the H~ chopping cycle and the resulting ^i^y was recorded using a slow 
current meter. As a result, for that work / was not measured during data collection but 
under simulated data collection conditions and an average value was used. Additionally this 
factor was used to extract I^l and I°^L from (I^° p ). 

For the present work we are also using a new calibrated neutral detector in combination 
with a fast current amplifier to record the H particle current at each phase in the chopping 
pattern and to monitor it throughout each data run. This modification is described in 
6j. Thus, during a data run, we are now able to directly measure Jr- which is needed in 
Eq. (j3J) . In our previous work, the H particle current was also monitored with a fast current 
amplifier; however, the neutral detector was not designed for absolute measurements. So 
to analyze those results, using the new detector we measured the H particle current due to 
photodetachment (PD), J™, to determine the neutral-to-anion (nta) ratio 

jPD 

/nta = j q q (6) 

under simulated data collection conditions. This factor, combined with the extracted I^L 



discussed above, was then used in [6| to determine /h for Eq. (jl]). 

We have also installed a BPM immediately before the neutral detector, at a distance 
of 2055 mm from the first BPM in the interaction region. Turning off the voltage of the 
LCD, allows the H~ beam to pass through the hole in the outer plate of the LCD. We used 
this additional BPM to measure the position of both the H and H~ beams and verified the 
alignment of the beams over a much longer lever arm than was previously possible using 
only the two BPMs in the interaction region. We find that the full angle between the beam 
axes measured here is in good agreement with that reported in 
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C. Pressure corrections 



Any H^~ formed in the gas cell can be destroyed by subsequent collisions with He in 
either the gas cell or the analyzer region. The resulting products are not transmitted by 
the electrostatic deflectors into the CEM, thereby reducing the apparent signal and rate 
coefficient. This small systematic shift to our data was overlooked in our previous work 
3, (]]. Here we quantified this minor correction for both our previous and present results. 

We measured the H 2 attenuation using an approach similar to the one we used to de- 
termine the He gas cell column density in 

as 

. Reconfiguring the ion source to produce 
H^~ and the apparatus to transmit H 2 beams, we used the electrostatic quadrupole after the 
interaction region to direct the beam into a Faraday cup where we measured the unatten- 
uated H^~ current, + . We then guided the beam through the gas cell and measured the 
transmitted current, J H +, on the outer plate of the UCD. With no He in the gas cell, the 
UCD reading was over 95% of that in the Faraday cup. The measured attenuated data were 
corrected for this slight difference in the unattenuated current readings. 
The H^" attenuation as a function of gas density is given by 

exp(-cr d iVHe), (7) 



1° 

where <7d is the total H^~ destruction cross section and Nu e is the helium column density. 
Following the methodology of [6], the column density can be expressed as 

N Re = / 7lHe(0*+ / 7lHe(0*+ / n He (l)dl. (8) 

J quad J gas cell J analyzer 

Here nn e (0 is the helium density and dl the infinitesimal path length. Using the same model 
as [(J, we take the pressure to be constant in each of these regions and re-express Eq. [H]as 

Nr c = n t li + n 2 l 2 + n 3 l 3 . (9) 

The He density in the quadrupole is n\ and the path length l\ = 5.0 ±1.0 cm. In the gas cell 
the He density is n 2 and the path length l 2 = 78.7 ± 1.0 cm. The He density in the analyzer 
region is n 3 and the path length l 3 = 35.4 ± 1.0 cm is the distance that the ions travel before 
striking the UCD. All uncertainties here and throughout the paper are given at an estimated 
la statistical confidence level. The respective densities were calculated from the measured 
pressures using the ideal gas law at the laboratory temperature which was stabilized at 293 K 



for both the work of and our new results here. The ratio of the measured pressures 

in each section were P1/P2 = 0.137 ± 0.019 and P3/P2 = 0.105 ± 0.034. The uncertainties in 
these ratios are due to the manufacturer-quoted accuracies of the pressure gauges (10% for 
Pi and p 2 and 30% for p 3 ). 

Attenuation data were collected for pressures up to ~ 4.5 x 10~ 4 Torr and are shown 
in Fig. [TJ From a fit to these data we extracted a cross section of (2.75 ± 0.29) x 10~ 16 
cm 2 at an energy of 10 keV amu -1 . This estimated uncertainty is due to the error in 
the attenuated and unattenuated current readings (3% each) and the uncertainty in the 
He column density (10%). The error in this latter quantity was estimated by adding the 
uncertainties from each segment nj/j of the total column density. The errors in the path 
lengths and gas densities (i.e., pressures) have been given above. 

Collisional destruction of has also been studied by [13] who reported cross sections of 
various outgoing channels for ion energies from 2-8 keV amu -1 . We have derived a total 



destruction cross section by summing the relevant channels in [13j. Those results, shown in 
Fig. El indicate that the cross section is essentially constant between 2 and 8 keV amu -1 . 
Our result at 10 keV amu -1 , also shown in Fig. El is in good agreement with this trend. 

To determine the expected signal attenuation factor and correct for the signal loss 
we use our measured Hj~ destruction cross section combined with Eq. ([7]). The appropriate 
He column density is given by 

K e =~n 2 k + n 3 l' 3 , (10) 

where the factor of 1/2 takes into account that on average the H^~ ions will be formed in 
the center of the gas cell and 1' 3 = 57.9 ±1.0 cm is the distance from the end of the gas cell 
to the CEM mouth. Using these values we calculate from Eq. [7| that the signal attenuation 
with 2 x 10 -4 Torr He in the gas cell is 0.92 ± 0.01. The signal must be divided by this 
factor to correct for the attenuation. This corresponds to an (8.6 ± 1.2)% upward shift in 
the data. The uncertainty in this correction is estimated by propagating through Eq. [7] the 
quadrature sum of the uncertainties from both N{± e in Eq. ( flOl) and cr^. 

III. UNCERTAINTIES 

The various systematic uncertainties for the measurement are given in Table [B Values 
are listed at an estimated lex statistical confidence level. We have grouped them into two 
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sets. The errors listed in the top half of the table add in quadrature to ±12% for each 
data point. This represents the relative uncertainty between our old and new data sets and 
also at different energies within each set. Adding this in quadrature with the remaining 
uncertainties in the bottom half of the table yields the total systematic error of ±24%. A 
detailed discussion is given in {3, (| for the various uncertainties not already discussed here. 



IV. THEORY 



A. Earlier calculations 



In our previous work, the AD cross section was calculated using non-local resonance 
theory and considering only the coupling of the H ± H~ and H2 + e~ channels through 
the lowest metastable state of 2 £+ symmetry (see 3), Q] for details). This state is one 
of two connected to the H + H asymptote (not counting the spin degeneracy). Potential 
energy curves for both states are shown in Fig. [3J The second state of the 2 E+ symmetry 
is repulsive and usually neglected in the calculations. The validity of this approximation 



non-local calculations 



is supported by the very good agreement between our experimental results [3_[ [f| and our 



below 1 eV, even after the ~ 9% pressure correction of the FT 



signal described in Sec. Ill CI which was not accounted for in [3|, [6j. 



B. New calculations 



We have extended our experimental results to ~ 5 eV, entering a regime where AD 
via the 2 £+ state becomes possible. Figure [3] shows that for sufficiently large energies the 
colliding H + can penetrate into the autodetachment region along the repulsive 2 £g 
state. This region is defined as the range of internuclear separations R where an electron 
can escape the anionic system, i.e., the potential energy curves of the H~~ system are above 
those for neutral H 2 . This occurs for the 2 S+ state at R < 5 do, where ao is the Bohr radius. 
Particles colliding along this state can penetrate into the autodetachment region for energies 
> 0.75 eV motivating calculations for AD via this state. 

Due to the different symmetry of the molecular orbitals, the 2 E+ and 2 Eg contributions 
to the AD cross section can be calculated separately. Thus we need only carry out new 
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calculations for the 2 S+ state. A brief description of our approach is presented below, using 
atomic units. A more detailed discussion will be given in a future publication. 

n 

Non-local resonance theory is explained in detail by [20(. The main idea is as follows. The 
electronic state <pd, describing the colliding partners in the H + H~ channel, is diabatically 
prolonged to small R. It is also assumed to be coupled to the H 2 + e~ electronic continum 
states (pk through the matrix element 

Vdk(R) = {4>d\H e i\(j>k) 

where H e \ is the electronic Hamiltonian. The non-local resonance model is parametrized 
by three functions: Vo(-R), Vd(R), and Vdk(R)- The potential energy curve for the neutral 
molecule Vo(-R) and for the anion Vd(R), are functions only of R. The coupling element 
Vdk{R), however, depends on both R and the momentum of the detached electron k. 

Once Vo(-R), Vd(R), and Vdk(R) are known, the electronic dynamics of the system is 
fully parametrized and the nuclear dynamics can be treated as a motion in the non-local 
energy-dependent effective potential 

V d (R) + J V dk (E - \k 2 -T N - V (R) + ie)- l V* k kdkdn k . 

T/v is the kinetic energy operator for the nuclei, dQk is the differential solid angle for the 
outgoing electron, and e is the usual positive infinitesimal of scattering theory. We solve the 
nuclear dynamics and calculate the cross sections using the method of Q]. 

In order to include <pd for the 2 S+ state, we have to fix the parameters of the non- 
local resonance model for this state. The proper procedure for calculating these parameters 
involves extracting the discrete state 4>d from the continuum <pk, employing the projection- 



operator technique. This procedure was followed in [21] for the 2 E+ state and we used it as 
an input for our calculation But it is also possible to fix the model parameters by fitting 
the fixed nuclei scattering data. We follow this latter procedure here. 

To fix the coupling amplitude, we assume the separable form Vdk(R) = g(R)f(k), where 



the k dependence is determined by the Wigner threshold law 22] with an exponential cut-off 

f(k) ~ k 2l+1 e- ak \ (11) 

where a is the cut-off parameter. The angular momentum I value in Eq. (fTTT) is given by the 
lowest electron partial wave allowed by symmetry (discussed below). The R dependence is 



determined from the calculated local decay widths T(R) of [15]. The potential energy curve 
'or the anion Vd(R) is constructed from 



,151,1161 and extended to larger R using the data of 
~y\ . The data for the potential near the crossing of the neutral and anion potential energy 
curves are missing. Nevertheless the analytic behavior near the crossing has been discussed 



in detail by 20j. With this knowledge, the potential energy curve can be interpolated 



through the crossing as has been done before for hydrogen halides [231 ] . The actual shape 
of the Vd(R) and Vq(R) crossing is modified by the interaction of the electron scattering 
continuum with the threshold behavior given by the Eq. (fTTj) . Vq is from 18 1. 

The decay of the odd symmetry anion 2 £+ state to the even neutral 1 El_state is possible 
only through release of an electron with odd angular momentum. In 0,1^] we considered 



only I = 1 (p-wave scattering) since the calculations of 24j show that the next allowed 
I = 3 contribution is suppressed by almost two orders of magnitude for the energy range 
of interest. For the anion 2 £+ state decaying to the neutral X S+ state, the symmetry 
remains unchanged, requiring release of an electron with even angular momentum. Here we 
considered only I = (s-wave scattering). In each CclSC, clS Eq. [U] shows, Vdk is strongly 
suppressed for higher angular momenta at the k < 1 values in our experimental results. 

The anion 2 £g state can also decay to the first excited 3 E+ state of H 2 . These states are 
both repulsive and lie very close together. This decay, however, requires an odd value for 
I. With I = 1 and k < 1, f(k), and hence Vdk, is strongly suppressed compared to the 2 S+ 
to decay channel with I = 0. The effect from the transition between repulsive states is 
thus expected to be small at low energies and was not included here. 

Once the model parameters are fixed, <tad can be calculated using the methods described 



in jg]. Figure S] shows our results. As expected from the previous good agreement of our 
experimental and theoretical results, the new contribution is small and notable only for 
E r > 0.75 eV. This is the threshold where the colliding particles overcome the barrier in the 
repulsive interaction potential and penetrate into the autodetachment region. 

Both the 2 S+ and 2 S+ contributions decay rapidly to zero for energies above ~ 1 eV. 
This is due to the competing process of collisional detachment 

IT + H^H + H + e-, (12) 

which opens up for E T = 0.76 eV and wins at higher energies. This is discussed in Sec. I VI I 
from the point of view of general energy conservation arguments. 
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C. Contributions from quasi-bound H2 states 

At high angular momentum (J > 10), the colliding H~ and H systems can autodetach 
into quasi-bound H 2 . These states, sometimes referred to as orbiting or shape resonances, 
lie above the separated atoms limit for H + H. Such high J levels, temporarily stabilized by 
the centrifugal barrier, will eventually dissociate spontaneously and are therefore generally 
not considered in AD cross section calculations. However, the lifetime for all but a few of 
these resonances well exceeds the flight time from the interaction region to the gas cell and 
so most are expected to contribute to the experimental signal. 

The H 2 flight time from the interaction region to the gas cell is (737 ±640) ns. The mean 
is the center-to-center distance, the upper limit is from the start of the interaction region 
to the end of the gas cell, and the lower limit from the end of the interaction region to the 
start of the gas cell. Quasi-bound H 2 (i.e., in high J levels) can strip in the He gas cell and 
will form H^ in similarly high J levels. As the H^ potential supports stable ro- vibrational 
levels up to J = 35, we assume that any such H^ formed will be stable and will reach the 
detector. 

In order to compare to our measured results, we have added the contribution of these 
quasi-bound H 2 states to our calculations for AD via the 2 E+ state. So as to mimic the range 
of experimental lifetimes, we have investigated the effect of cutting out states with lifetimes 
less than 100, 700, and 1400 ns and found no significant differences. In the end we included 
contributions from all resonances with lifetimes longer than 700 ns. The contribution of these 
states is of comparable size to the 2 £+ state contribution. The effect of these resonances 
for AD via the 2 S+ state has not been considered as that would be a small correction to 
an already small contribution. Lastly, we note that the significance of these resonances for 
molecular hydrogen formation in plasma environments will depend on whether the states 
can relax to stable states of H 2 before they dissociate by tunnelling. 



V. RESULTS 



Relative energies E r are controlled by varying the potential of the floating cell U{. In 
{3, [s] data were collected for E r < leV (\U{\ < 281 V). Here we have extended the energy 
range to E T < 4.83 eV (| Uf | < 621 V). Data are collected by stepping U{ in voltage. The 
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present work uses voltage ranges smaller than our earlier measurements. For \U{\ < 441 V, 
Uf was scanned across 60 V ranges in 10 V steps and for \Uf\ > 441 V the scanning was 
across 120 V ranges in 20 V steps. 

Our measured rate coefficients for reaction ([T]) are plotted in Fig. as a function of 
average collision energy (E T ) < 1.0 eV. The black circles represent our new results and 
the red triangles our previous work. Both have been corrected for the attenuation of the 
ions. The error bars on each data point display the la statistical uncertainty. There is 
an additional ±12% relative systematic error on each data point which is not shown. The 
good agreement between our new and previous results indicates that there were no hidden 
systematic errors due to our previous inability to measure and monitor / and / nta during 
data acquisition. 

A final potential source of systematic error which we investigated was to verify the linear- 
ity of the gas stripping method used to convert the product H 2 molecules into the measured 
signal. Here we measured the AD rate coefficient as a function of helium gas cell pres- 
sure for (1 — 3) x 10~ 4 Torr. Table [TT] shows the results of these AD measurements at 
(E T ) = 16 meV versus pressure. Taking into account the attenuation of the signal ions, 
to within the uncertainties the data show no dependence on gas cell pressure. 

Given the good agreement between our results in {3, 6] and our new data, we have 
merged them together using a statistically-weighted averaging method. We also included 
our pressure test results in this average. The la counting statistics of each data point 
were used for the weighting. All data sets were also measured on the same relative energy 
grid. Figures E] and [7J show the averaged data for (E T ) < 1 eV plus the new data we have 
collected for 1.0 eV < (E T ) < 4.83 eV. Also shown are the cross section calculations of {sUsl, 
supplemented by our new theoretical work here, multiplied by v T , and convolved with the 
experimental energy spread. Figure [7J shows the theoretical results with and without the 
effects of the H 2 orbiting resonances included. As is clear from the figures, we find good 
agreement with theory throughout the measured energy range. The contribution due to 
orbiting resonances of H 2 can also be seen in Fig. [7J as the experimental data are shifted to 
slightly higher energy compared to calculations which do not include these resonances. 
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VI. DISCUSSION 



The good agreement that we find here both with our previous results and with our up- 
dated theory strengthens our confidence that theory and experiment have finally converged 
for reaction ([T]). Including AD via the repulsive 2 S+ state increases the cross section by 
an amount smaller than we are currently able to measure experimentally. The resulting 
theoretical thermal rate coefficient is only 1.3% larger than that for only the attractive state 
at temperatures of 4,000 K, 3.5% at 8,000 K and 4.4% larger at 10,000 K. These are sig- 



nificantly smaller than the ~ 25% experimental accuracy with which we have been a 
benchmark theory. Hence, we continue to recommend the thermal rate coefficient of 
modeling plasma temperatures below 10 4 K. 



e to 
$ for 



Additionally, our results continue to imply that the reason for the discrepancy seen with 
the results of {4] lies in the data of {t| used for normalization. This is further supported by the 



theoretical AD work on hydrogen halides of 



25j . They used the same theoretical approach 



as we do here and found systematically higher AD rate coefficients than the experimental 
work of Q] . It appears to us that a re- measurement of reaction (J2J) using a technique different 
from that of j3| is clearly called for to resolve this dilemma. 



Our results also verify the predictions of [8j and our new work here that the AD cross 
section for reaction ([I]) should decrease to essentially insignificant values for E T > 2 eV, as 
shown in Fig. |6j A simplified adiabatic description of the AD reaction can provide good 
insight into the physics behind this prediction. We consider here only the 2 £+ symmetry. 
Similar arguments can also be given for the 2 S+ state. 

Initially the H~ and H approach one another along the attractive 2 £+ electronic state. 
This state crosses into the autodetachment region at R ~ 3 ao- Adiabatic theory dictates 
that the system remains electronically in the ground state. Inside the autodetachment region 
the ground state is the x £g state of neutral H 2 plus a free electron with zero kinetic energy. 
Conservation of energy requires that the final state energy equals that initially available 

E v = E T + D - E EA . (13) 

Here E v is the excitation energy of the vibrational level v formed in the process; D is the 
4.48 eV dissociation energy gained by formation of H 2 in the v — vibrational and J = 



rotational level 26]; and Eea = 0.76 eV is the electron affinity required to neutralize the H 
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and form H (19J. For E r > Eea, the system lies in the dissociation continuum (E v > D ), 
resulting in the formation of H + H + and not H2 + e~. 

In reality the AD process is not exactly adiabatic. This is manifested by the release of 
an electron with a nonzero kinetic energy E e and we can rewrite ([TBI as 

E v + E c = E r + D - E EA . (14) 

The nonadiabatic exchange of energy between the electron and protons is weak though; 
detached electrons do not have a large kinetic energy. Our full calculations for the 2 E+ state 
8J show that only a negligible amount of electrons can have energy above ~ 1.5 eV. Taking 
into account that the largest possible value of E T will occur for E v = Do, this leads to the 
prediction that the AD process will cease for E r > Eea + 1-5 eV. For reaction ([1]), this 
corresponds to E t > 2.26 eV. A similar argument has been suggested for the decrease in the 
cross section for protonium formation in collisions of antiprotons with hydrogen atoms (see 
[27I ] for a review). Note that we have ignored the insignificant kinetic energy of the final H2 
molecule E^ 2 , as conservation of momentum gives E H2 = (m c /m H2 )-E e <C E c . 

Continuing the protonium analogy, one would expect a sharp decrease in a ad immediately 
after the collisional detachment threshold at E x = 0.76 eV. In the H~ + H collision, the 
drop in the cross section occurs at higher energies. This is related to the threshold law given 
by Eq. ( TTTT) with / = 1 for the dominant ungerade channel. As a result, the coupling Vdk 
vanishes for zero detached electron energy and rises smoothly as the energy increases. The 
electron energy in Eq. ffT4"|) thus can not be exactly zero, but remains relatively small. The 
smooth decrease in a ad above 1 eV, confirmed by the present experiment, thus provides a 
good test of the theoretical description of the electron release amplitude. 

Lastly, the decrease of the AD cross section is indeed slightly weakened by positive con- 
tributions of the 2 S+ state as shown in Fig. H] and orbiting resonances as seen in Fig. [71 
However, the decreasing trend above leV, controlled by Vdk, is still dominant (e.g., Fig. [6]). 



VII. SUMMARY 

We have modified the experimental methods used in j^, 5, 6] 
to E r < 4.83 eV. Additionally, we have performed several modifications to better control 
potential systematic errors. We find good agreement between our previous and new data 



6] to measure reaction (JTJ) up 
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sets. To within the experimental uncertainties, we also continue to find good agreement 
with the calculations of {3, 8] which have been extended here to include contributions from 
the repulsive 2 S+ H^" state and for the effects on the experimental results due to orbiting 
resonances of H 2 for E T > 0.76 eV. In particular, we confirm the predictions of [8] that this 
reaction turns off for E r > 2 eV. Similar behavior has been predicted for the formation of 
protonium from collisions of antiprotons and hydrogen atoms 27]. 
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1 2 3 4 5 

He pressure (10~ 4 Torr) 



FIG. 1: Attenuation of the ion beam as a function of helium gas cell pressure. The circles 
represent the statistically-weighted mean from three sets of measurements. The error bars are 
smaller than the plotted circles. The line shows the best exponential fit. 
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Collision energy (keV amu ) 



FIG. 2: Experimental cross sections for total destruction versus ion beam energy for + He. 
The open squares are the results of [13] while the circle represents our measurement. The error 
bars for each data set give the total la experimental uncertainty. 
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FIG. 3: H 2 and H2 potential curves versus internuclear distance in units of the Bohr radius qq. 



The H 2 attractive 2 S+ electronic state is given by the dashed curve 



electronic state by the dotted curve constructed using the data of 



of [l7] above ~ 5 «o. The separated atoms limit (SAL) ibr these two pot rutin] energy curves is 
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14] and the repulsive 2 £g 
16[ below ~ 4 ao and those 



18(] with an SAL of 



H ( S) + H( S). The solid curve shows the H2 Sg electronic state from 
H( 2 S) + H( 2 S). The energy difference between the two limits is determined by the H electron 
affinity energy E EA = 0.76 eV Q- 
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FIG. 4: Theoretical cross section for H~ + H — > H2 + e~ as a function of the relative collision 
energy E r . The dashed curve shows the results via the attractive 2 S+ state, the dotted curve 
via the repulsive 2 £+ state, and the solid curve the sum of the two. 
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(E r ) (eV) 

FIG. 5: (Color) Experimental rate coefficient (dAD^r) as a function of the collision energy (E T ). 
The black circles show our new results and the red triangles our previous results of Q, [f| corrected 
for the H^" attenuation. Although our new results extend up to (E T ) < 4.83 eV, here we show 
only up to the maximum {E T ) of our previous results for comparison. The error bars show the la 
statistical uncertainties. There is an additional ±12% relative systematic error on each data point 
which is not shown. 
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FIG. 6: The circles show the statistically-weighted mean of the experimental rate coefficients 
(cjAD^r) from our previous jj, 6] and current work as a function of the collision energy (E r ) (see 
text). For (E T ) > 1.0 eV the data are solely from the current measurement. The error bars 
represent the la statistical uncertainties. The solid line is from the cross section calculations of 
3, 8], supplemented by our new theoretical work here, multiplied by v r , and convolved with our 
experimental energy spread. The effects of the H2 orbiting resonances have been included in the 
calculations shown here. 
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(E r ) (eV) 

FIG. 7: Same as Fig. [6] but on a linear scale. The dotted curve shows the calculations without the 
effects of the H2 orbiting resonances included. 
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TABLE I: Summary of systematic uncertainties at an estimated la confidence level. Uncertainties 
are treated as random sign errors and added in quadrature. 

Source Error (%) 

Background subtraction 5 

Anion current 3 

Neutral current 10 

Beam overlap 3 

Total relative errors from above 12 

Stripping cross section 16 

Effects of unknown rovibrational population 10 

Analyzer transmittance 1 

Grid transmittance 1 

CEM detection efficiency 2 

Overlap length 1 

Helium gas cell column density 7 

Attenuation 1 

Total systematic uncertainty 24 
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TABLE II: Rate coefficient results at (E T ) = 16 meV versus helium gas cell pressure. Our theo- 
retical results are also shown for comparison. 



Pressure (40 ^ Torr^ 




Rate Coefficient (10 9 cm 3 


s- 1 ) 


Value 


Statistical Uncertainty 


Relative Uncertainty 


1.0 


5.6 


± 0.5 


± 0.7 


2.0 


5.2 


± 0.4 


± 0.6 


3.0 


5.2 


± 0.4 


± 0.6 


Theory 


5.0 
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